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ABSTRACT 

Near-IR and X-ray flares have been detected from the supermassive black hole Sgr A* at the center of our 
Galaxy with a (quasi)-period of ~ 17-20 minutes, suggesting an emission region only a few Schwarzschild 
radii above the event horizon. The latest X-ray flare, detected with XMM-Newton, is notable for its detailed 
lightcurve, yielding not only the highest quality period thus far, but also important structure reflecting the 
geometry of the emitting region. Recent MHD simulations of Sgr A*'s disk have demonstrated the growth 
of a Rossby wave instability, that enhances the accretion rate for several hours, possibly accounting for the 
observed flares. In this Letter, we carry out ray-tracing calculations in a Schwarzschild metric to determine 
as accurately as possible the lightcurve produced by general relativistic effects during such a disruption. We 
find that the Rossby wave induced spiral pattern in the disk is an excellent fit to the data, implying a disk 
inclination angle of ^ 11° . Note, however, that if this association is correct, the observed period is not due 
to the underlying Keplerian motion but, rather, to the pattern speed. The favorable comparison between the 
observed and simulated lightcurves provides important additional evidence that the flares are produced in Sgr 
A*'s inner disk. 

Subject headings: accretion — ^black hole — Galaxy: center — instabilities — magnetohydrodynamics — ^relativity 



1. INTRODUCTION 

Long-term monitoring of Sgr A* has thus far res ulted in 
the discov ery of several near-IR dGenzel et alj|2003l) and X- 
ray flares (iBaganoff et alJ l2001h from this object. One of 
the most recent events, observed with XMM-Newton in 2004, 
has produced compelling evidence of a significant modulation 
in the X-ray lightcur ve with a quasi-period ~ 22.2 minutes 
(iBelanger et alJl2006r) . Similar periodic fluctuations had also 
been seen earlier in the near-IR events. 

A period of ~ 22.2 minutes is rather intriguing because 
simple considerations dMelial 2001) would place the corre- 
sponding emission region at roughly 3 Schwarzschild radii 
(r, = IGMIc^) fo r a black hole (BH) mass of ~ 3.4 x lO*" 
Mq (ISchodel et al.ll2003i) . It is not yet understood, however, 
why the Keplerian radius corresponding to this period is ac- 
tually less than 3rs. It could be that the difference is due to 
the spin of the BH, which moves the inner disk radius in- 
wards for prograde rotation. Or it could be that the so-called 
stress edge — the location where the inspiralling material ac- 
tually detaches from the rest of the magnetized disk — is dy- 
namically impo rtant in es tablishing where the disk e mission 
terminates (.Krolik & Hawle y 2002: lMelia et al.ll2007h . 

Part of the uncertainty is due to the fact that the struc- 
ture of the disk surrounding Sgr A* is itself not fully un- 
derstood yet. We have recently begun to simulate the be- 
havior of the hot, magnetized disk during a disruptive event 
that may be responsible for the flares, under the assump- 
tion that the instability is induced by low angular momen- 
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tum clumps of plasma "raining" inwards towards the Kep- 
lerian region (Tagger & Melia 2006; see also Chan et al. 
2006). This scenario is motivated by extensiv e hydrody- 
namic (iMeUa & Cokei'll 999t ICuadra et al.ll2005h and MHD 
(llgume nshchev & NarayanI l2002h simulations which show 
that for the stellar-wind fed conditions at the Galactic center, 
the average specific angular momentum of gas captured grav- 
itationally by Sgr A* is too small to sustain a 'conventional' 
(i.e., typically large ~ 10^' r^) disk. Instead, only clumps 
of plasma with relatively small angular momentum venture 
inwards and merge with — essentially, 'rain' onto — the com- 
pact disk at the circularization radius , which for Sgr A* is 
< 10 - 10^ r, (see also Melia & Falcke 2001; Melia 2007). 

These MHD simulations (Tag ger & Melia.2006 ) did indeed 
establish the result that the merger of an infalling blob with 
the existing Keplerian flow induces a Rossby wave instability 
(RWI) leading to the total disruption of the disk on a timescale 
relevant to the Sgr A* flare phenomenon. But unless the en- 
suing X-ray modulation is due to a QPO response in the in- 
spiralling matter (see Chan et al. 2006), one cannot escape 
the fact that general-relativistic effects are essential in produc- 
ing a periodicity in the lightcurve due to radiation by an az- 
imuthally asymmetric emitter. The pattern of modulation over 
one complete cycle would be the result of several influences, 
including a Doppler shift, light- bending, and lensing e ffect 
near the BH's event horizon (e.g.. lHollywood et al.lll995T) . 

In this Letter, we ta ke the results from the (non-relativistic) 
MHD simulation of Tagger & Melia (2006) and carry out a 
full ray-tracing calculation of the lightcurve produced by the 
disruption for direct comparison with the XMM-Newton data. 

2. OBSERVED LIGHTCURVE OF AN X-RAY FLARE FROM SGR A* 

An inspection of Fig. 5 in lBelanger et al.l (l2006l) shows that 
the effects of gravitational light-bending, lensing, Doppler ef- 
fect and travel time delay may have helped to shape the folded 
lightcurves from the 2004 August 3 1 event. One cannot yet 
discount the possibility that an actual periodic dynamical ef- 
fect may have also contributed to the modulation seen in the 
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emissivity, but this is necessarily model-dependent and the 
signature may not be unique. For example. Tagger & Melia 
(2006) found that a quasi-periodic modulation could be ex- 
cited by non-linearities in the evolution of the spiral-Rossby 
pattern, which would sit on top of the modulation one would 
see due to general relativistic effects. For the purposes of this 
Letter, we will adopt the simplest assumption — that the mod- 
ulation is due predominantly to general relativistic effects. 

The 2-10 keV folded curve from the 2004 August 31 event 
(see Fig. O looks broader and, except for one significant da- 
tum standing well away from the rest, looks more like a con- 
tinuous, relatively smooth modulation rather than the sharp 
changing profile produ ced by an orbiting hot spot (see e.g., 
iHollvwood et al.iri995h . It is partially for this (phenomeno- 
logical) reason that the Rossby wave instability is promising 
and worth investigating further here, because the disruption it 
causes is global, offering the possibility of producing a grad- 
ual modulation in the lightcurve when general relativistic ef- 
fects are included. 

3. NUMERICAL SIMULATION OF THE INSTABILITY 

The instability we have simu lated with our MHD code ha s 
a long history, dating back to iLovelace & Hohlfeldl (Il978l) . 
who showed that a disk presenting an extremum of a quantity 
X (later dubbed vortensity) was subject to a local instability 
of Rossby vortices. The requirement of an extremum is sim- 
ilar to that giving rise to the Kelvin-Helmho l tz in stability of 
sheared flows. More recently, Lovelace et al.l (11999 ) renamed 
it the Rossby Wave Instability (RWI) and developed the the- 
ory and numerical simulation. 

In isothermal, unmagnetized disks, X is the specific vortic- 
ity averaged across the disk thickness, 

vxy ^ ^ 

where S is the disk's surface density, Q its rotation frequency, 
and - 4Q^ -i- 2Q.Q.'r is the epicyclic frequency squared. 
The extremum of X appears to be due to an extremum in 
the radial density profile. To understand how the instability 
is driven, we note that Rossby waves in disks form normal 
modes trapped near the extremum of In the MHD form 
of the RWI, the disk is threaded by a vertical (poloidal) mag- 
netic field BQif). Its properties are essentially the same as 
those discussed above, except that here the critical quantity 
is JLb - K^^/imBf.^), and the growth rate can be higher be- 
cause of the long-range action of the Lorentz force coupling 
the Rossby vortices. 

4. RAY TRACING CALCULATIONS 

A typical profile of the inner disk during the Rossby wave 
growth is shown in Fig. 3 of Tagger & Melia (2006). In this 
paper, we present the lightcurve and images associated with 
this disruption, from the vantage point of an observer at infin- 
ity. The calculation is carried out with a fully general relativis- 
tic ray-tracing code. The RWI arises in the compact accretion 
disk surrounding the (Schwarzschild) BH, and we describe 
its morphology using coordinates in the co-rotating frame 
(r, 6, if). The modeled accretion disk is thin and the RWI may 
be considered to lie in the equatorial plane {9 - n/2) of the 
compact object. The observer is located at infinity with view- 
ing angle / relative to the z'-axis in the non-rotating frame, at 
(observer) polar coordinates (r', 6', (p'). The deflection angle 
of a photon emitted by plasma in the Rossby-unstable region 
is ijj, and varies periodically with cos i/r - cos / cos if. 



These emitted photons are deflected by the BH and intersect 
the observer's detector plane at infinity. The distance between 
the line-of-sight and the point at which the photon reaches the 
detector is defined as the impact parameter b. Using this ge- 
ometry, the deflection angle of the photon's trajectory may be 
obtained with the light-bending relation between a (the angle 
between the emission direction of the photon and the direc- 
tion from the center of the BH to the location of the emitter) 
and ijj from the geodesic equation 

*^rfl^-;^('-7)l- « 

This procedure yields the impact parameter b - r{\ - 
rslrY^ sin a of the photons in terms of the emitting radius 
r, and ultimately allows us to calculate the flux at infinity. 
A detailed desc ripti on of this geometry i s provided in, e.g., 
iLumined (Il979l) . and lFalanga et al.l (l2007h . In our derivation, 
we shall use the same notati on and geometry described in 
iPoutanen & Gierliriskil ( l2003h . though pertaining to a local- 
ized region (or clump) in the disk. The system of units is cho- 
sen such that G - c - 1 ; in this standard coordinate system, 
the BH's horizon occurs at the Schwarzschild radius - 2M, 
where M is the mass of the compact object. 

In the simulations we report here, the emitting region is ge- 
ometrically thin, and generally optically thin as well. Our rays 
include an integration through the whole emitting depth, but 
because most of the emissivity is concentrated near the disk's 
plane, we effectively have a situation in which the rays them- 
selves appear to begin very close to this plane. Rays leaving 
the disk in directions that eventually take them around the BH, 
heading toward the observer, contribute much less to the over- 
all flux and we ignore them here. 

The general-relativistic effects to be considered are now: 
(/) light-bending (see above), (//) gravitational Doppler effect 
defined as (1+z), (///) gravitational lensing effect, t/Qobs - 
bdbdiflD^ (with D the distance to the source), expressed 
through the impact parameter, and iiv) the travel time delay. 
We calculate the relative time delay between photons arriv- 
ing at the observer from different parts of the disk, using the 
geodesic equation. The first photon reaching the observer 
is the photon emitted at phase = 0, on the closest radius 
(f — rmax) to the observer We set this reference time, Tq, 
equal to zero. 

The observed flux at energy E' is FobsiE') = Iohs(E')dnobs, 
where Iobs(E') is the radiation intensity observed at infin- 
ity and t/Qobs is the solid angle on the observer's sky in- 
cluding relativistic effects. Using the relation Iobs{E',a') - 
(1 +zY^hm{E, a), a Lorentz invariant quantity that is constant 
along null geodesies in vacuum, the intensity of a light source 
integrated over its effective energy range is proportional to the 
fourth power of the redshift factor, /obs(Q'') - i^+zy^hmir^ip), 
hmif, ^ ) being the intensity measure d in the rest frame of the 
clump ( iMisner. Thorne. & Wheeledll973h . 

Although our magnetohydrodynamic simulation, and corre- 
sponding ray-tracing calculation, are quite sophisticated, they 
are nonetheless still somewhat restricted in that we have not 
allowed for a completely self-consistent coupling between the 
plasma and the radiation. The infalling plasma radiates inef- 
ficiently, so this is not a serious deficiency as far as the dy- 
namics is concerned. Furthermore, properly modelling the 
compression of the gas would require fully 3-D simulations, 
which MHD codes are still unable to handle in the conditions 
(a disk threaded by a near-equipartition vertical field) that we 
use. For simplicity, the MHD simulation was thus carried out 
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Fig. 1 . — Horizontal: Four snapshots of the disk (spiral pattern) as it would appear to an observer looking along a line-of-sight inclined by 30° (top), 60° (middle), and 80° (bottom), 
relative to the disk's symmetry axis, at phases if = 0, jijl, n, and 3;r/2, respectively. Variation in each column is due to general relativistic light-bending, lensing, and Doppler shifts, 
including the relativistic time delay for the various inclinations. The lightcurves attached to their respective images show how an observer views these overall eifects from infinity (see 
also Fig. 2). Slight differences between the two arins, resulting from the MHD siinulations, result in two slightly different peaks in the lightcurve during one full rotation of the pattern. 



assuming isothermal conditions. However the electron tem- 
perature, which dominates the radiative emission, responds 
to the compression of the gas and we model this here with 
an ad-hoc prescription, from the compression obtained in the 
simulation. This is sufficient to provide the main result of the 
present work, which is the form of the light curve that the 
spiral pattern and general relativistic effects generate. Also, 
we have not included non-local effects arising, e.g., from in- 
verse Compton scattering. In determining the surface emis- 
sivity, it therefore makes sense to take a simplified approach 
in which we include principally the parameter scalings, rather 
than their absolute values. This procedure will give us correct 
amplitudes in the lightcurve, though not the absolute value of 
the flux per se. 

Using the perfect gas law for an adiabatic flow (since the 
gas is radiatively inefficient), we can invoke a polytropic equa- 
tion of state with y - 5/3 and write the temperature as 
T oc p^l^ . This assumes further that the radiation pressure is 
negligible and T is not so high that the particles are strongly 
relativistic. Although the disk is not in full hydrostatic equi- 
librium, we can still argue that on average, we should have 
for a steady thin disk the gas density p - 'L/H(r,z), where S 
is the column density, and H is the disk height. These give 
p oc £3/4^-9/8 (usjjjg f j-Q denote the radius in the equatorial 
plane) and T oc £i/2^-3/4 synchrotron emissivity is there- 
fore js oc Bn„i oc ILiTp) oc lLp^^l^\ where the nonthermal par- 
ticle energy is roughly in equipartition with the thermal. We 
argue that the plasma is fully ionized and resistivity is mini- 



mal, so B is frozen into the gas, which means that B oc E. We 
therefore infer that 

Now the X-rays are produced via inverse Compton scatter- 
ing from the seed photon number flux. Thus, with Lseed °^ 
r' js, where js is the synchrotron emissivity in units of energy 
per unit volume per unit time, the soft photon flux scales as 
the emitted power divided by the characteristic area. That is, 
^seed 0^ j s / = fjs, which is going to be roughly the same 
scaling as the seed photon density, so Wseed fjs £^/'*r"^/^. 
The inverse Compton scattering emissivity is therefore j,f oc 
nntWseed 1)^/2^-11/4^ Thus, jx-ray ~ jic, and the surface in- 
tensity is /em oc j jx-rayds oc jx-ra\H, which gives finally 
/em 2'^/"^r"'^/^. The core physics in this expression is the 
surfac e density £(r, (^), which we take from the MHD simula- 
tion in lTagger & Melial ( l2006l) . 

In order to evaluate the flux at a given azimuthal angle (f 
and radius r, we first compute numerically (^(a), and then cal- 
culate the Dopper shift, lensing effects, and finally the flux 
Fobs as a function of the arrival time, including all the travel 
time delays described above. The simulation is carried out 
on a polar grid with - 256 and - 128 points, extend- 
ing from Tin = 3rs to 9Qrs. In this pass through the problem, 
in order to clearly separate the general relativistic modulation 
from the complex dynamical behaviour observed during the 
simulation, we have chosen to take a representative snapshot 
and rotate it at the frequency measured in the simulation. A 
full treatment of the whole simulation will be presented in 
iFalanga et ai] (l2007h . 
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Fig. 2. — Panels (a-c) show the lightcurve observed at infinity, resulting from general 
relativistic flux modulations associated with a Rossby-wave-disrupted disk. The period 
arises from the two-spiral-armed pattern speed, and is roughly 50 minutes for this par- 
ticular MHD simulation. Panel (a) is for a disk inclination of 30°, (b) is 60°, and (c) is 
B0° , in a one-to-one conespondence with the rows of Fig. 1 . Panels (d-f) show the same 
lightcurve, though now folded on half of a pattern-rotation period, i.e., 25 minutes. Each 
period is repeted once for clarity. 

5. RESULTS AND DISCUSSION 

Our principal results are presented graphically below, and 
may be understood with reference to the simulated images 
shown in Fig. [T] These are intensity maps projected onto 
the plane of th e sky , for the ray-traced perturbed disk in 
iTagger & Melial (12006) . at three different inclination angles: 
30° in the top row, 60° in the middle, and 80° for the bot- 
tom row. The four columns are snapshots taken at 4 (equally- 
spaced) phases of one complete pattern revolution. The gen- 
eral relativistic distortions depend strongly on inclination an- 
gle, which we employ in our search for the best fit to the mod- 
ulation in the X-ray lightcurve. 

This variation is demonstrated quantitatively in Fig. [2] 
which shows the lightcurves corresponding to the three incli- 
nations illustrated in Fig. [T] Note that at small inclination, we 
detect a gradual, broad modulation, whereas for the higher 
inclinations, we begin to see the effects of a two-spiral arm 
emitting region. Interestingly, the period associated with the 
pattern rotation in the Tagger &Melia (2006) simulation is 
about 50 minutes. Though this calculation was not optimized 
to fit the observed period, the fact that we see a bimodal mod- 
ulation from the two spiral arms suggests that one cycle in the 
data may in fact correspond to half a revolution of the pattern. 
In this figure, we therefore also show the calculated lightcurve 
folded over half a pattern period, corresponding to about 25 
minutes. 

The true test of relevance for this simulation lies in a direct 
comparison between the data and the inclination-dependent 
lightcurves. We emphasized in the introduction that the var- 
ious relativistic effects produce a unique profile, not easily 
confused with other periodic modulations, most of which tend 
to be sinusoidal. As we demonstrate in Fig. [3] the shape of 
the calculated lightcurve, particularly its amplitude, is quite 
sensitive to the inclination angle, which again, is most easily 
understood with reference to Fig. [T] For example, this figure 
includes 3 curves, corresponding to inclinations of 70°, 77°, 
and 80°. The middle curve produces the best fit (x^^j = 1.1), 
and it should be noted that the correspondence to the data is 
excellent, not only in terms of the modulation amplitude, but 



also for the shape of the lightcurve itself. Note that the simu- 
lated lightcurve corresponds to the bolometric flux emission. 
We have therefore compared the theoretical curve with the 
observed 2-10 keV lightcurve, rather than individually in dif- 
ferent energy bands, as presented in Belanger et al. (2006). 
There are several new ideas that we can take awav from 
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Fig. 3. — Lightcurve of the 20 04, August 31 flare folded with a phase of 1330 s, 
in the 2-10 keV energy band from lBelanger et all f2006). The best fit model is shown 
by the solid line using an inclination angle of 77°. The dasehed line represent upper 
and lower limits using 80° (bigger amplitude) and 70° (smaller amplitude) inclination 
angles, respectively. 



this work. First, if the periodic modulation seen in some 
flares from Sgr A* by both IR and X-ray instruments is real, 
it is not at all obvious that the periodicity is due to an un- 
derlying Keplerian period. One must be very careful, there- 
fore, in over-interpreting these periods in terms of a BH spin. 
Second, the shape of the X-ray lightcurve in this particular 
flare is too broad for it to correspond easily to a highly lo- 
calized hot spot on the disk, which instead would produce a 
more strongly peaked p rofile like that investigated earlier by 
[Hollywood et al.l ( Il995h . Instead, our work here argues for a 
more global disruption in the disk, at least for some events, 
like the 2004, August 31 flare observed with XMM-Newton. 
In this regard, the driving mechanism is likely to be an in- 
fall of clumps of plasma that merge with the existing compact 
disk and induce a Rossby-type of instability. Finally, the ob- 
served power density spectrum shows that the (quasi)-period 
is quite clean, without any evidence that the emission region 
is spread over a large range in Keplerian periods. This has 
been taken as some evidence that the disturbance must there- 
fore be highly localized, probably near the marginally stable 
orbit. But as we have shown here, the disturbance need not 
be that close to the event horizon in order to produce a mod- 
ulation with a period of only ~ 20 - 25 minutes. One may 
still get a narrow peak in the power density spectrum, as long 
as the modulation is due to a pattern rotation, rather than to 
motion along an orbit. But this will only work as long as the 
pattern has multiple components, such as the two spiral arms 
we have modeled in this paper. 
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sity of Arizona. FM is grateful for the hospitality of the APC 
in Paris, where most of this work was carried out. 
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